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Abstract. - The Lorenz 1963 dynamical system is known to reduce in the steady state to a 
one- dimensional motion of a classical particle subjected to viscous damping in a past history- 
dependent potential field. If the potential field is substituted by a periodic function of the 
position, the resulting system shows a rich dynamics where (standard) diffusive behaviours, 
ballistc motions and trapping take place by varying the model control parameters. This system 
permits to highlight the intimate relation between chaos and long-time deterministic diffusion. 



Introduction. Chaotic dynamical systems are known to exhibit typical random pro- 
cesses behaviour due to their strong sensitivity to initial conditions. Deterministic diffusion 
arises from the chaotic motion of systems whose dynamics is specified, and it should be dis- 
tinguished from noise induced diffusion where the evolution is governed by probabilistic laws. 
Diffusive (standard and anomalous) behaviours have been observed in periodic chaotic maps 
(see, e.g., refs. ^| and references therein) and in continuous time dynamical systems || |j. 
The analysis of deterministic diffusion is relevant for the study of non-equilibrium processes 
in statistical physics. The major aim is to understand the relationship between the deter- 
ministic microscopic dynamics of a system and its stochastic macroscopic description (think, 
for example, at the connection between Lyapunov exponents, Kolmogorov-Sinai entropy and 
macroscopic transport coefficients, firstly highlighted by Gaspard and Nicolis ||). 

In the present paper we introduce a model of deterministic diffusion in a one-dimensional 
lattice, suggested by a classical mechanics interpretation of the celebrated Lorenz 1963 sys- 
tem ||. The steady state chaotic dynamics of the Lorenz system can indeed be recasted 
as the one-dimensional motion of a classical particle subjected to viscous damping in a past 
history-dependent potential field (see refs. 0, ||, and ref. || for an earlier preliminar analysis). 

We shortly recall that the (scaled) Lorenz dynamical system is given by 
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with r > 1. For 1 < r < r c , where r c = a (a + b + 3)/(er — b — 1), three fixed points exist: 
(0,0,0) (unstable) and (±1,±1,1) (stable). For r > r c all fixed points are unstable, and the 
Lorenz system can exhibit either periodic or chaotic behaviour on a strange attractor set (see, 
e.g., ref. ]To| for a comprehensive exposition on the subject matter). 

In the steady state, the system (flh can be reduced to a one-dimensional integro-differential 
equation for the x coordinate ^ 

'x + T)x + (x 2 — l)x — — a[x 2 — l]px , (2) 

where a = (2a /b) -1,0= [26/ (r - 1)] 1/2 , 17 = (<r + l)/[(r - l)6/2] 1 / 2 , and the time is scaled 
by a factor [(r — 1)6/2] x / 2 with respect to the time coordinate in eqs. (|l|). The square brackets 
in eq. (||) indicate the exponentially vanishing memory, which is defined, for any suitable time 
function f(t), by 

poo 

— ks 



[f] k (t) = k dse- ks f(t~ S ). 
Jo 

According to eq. (g), the Lorenz system chaotic dynamics corresponds to a one dimensional 
motion in a constant-in-time quartic potential U(x) — (x 2 — l) 2 /4. In the presence of friction 
(rj 0) the motion is sustained by a time dependent memory term, which takes into account 
the system past evolution. 

Although eq. (|^) has been deduced from the Lorenz system (p]), it can be generalized to a 
wider class of equations showing similar dynamical properties |j. Indeed, it can be usefully 
recasted in the form 

x + rjx + {q(x) +a[q(x)]p}&(x) = (3) 

where the prime indicates the derivative with respect to x. Equation (|^) is obtained for 
$(x) = x 2 /2 and q(x) = x 2 — l. The generalized equation (|3|) can be regarded as the description 
of the motion of a unit mass particle subjected to a viscous force — r\x and interacting with 
a potential field <f>(x) through a dynamically varying "charge" qt(x) — q(x) + a[q(x)]/3. This 
charge depends both on the instantaneous particle position x(t) and on the past history 
{x(t — s) I < s < 00}. It is just the coupling of [<7(x)]/3 with the fixed potential field Q(x) 
the origin of an endogenous forcing term which can sustain the motion even in the presence 
of friction: the chaotic behaviour can actually arise from the synergy between this term and 
the viscosity. 

Moreover, one can easily verify that eq. ([|) corresponds to the generalized Lorenz system 

- a(y — x) 

-.-y + x+(r-l)(l-z)&(x) (4) 
= -bz + b[\q'(x)(y - x) + q(x) + 1] . 

The specific Lorenz model can thus be viewed as singled out from a quite general class of 
dynamical systems which can exhibit chaotic behaviour, their common essential property 
being an exponentially vanishing memory effect together with a viscous damping. 

In our previous works 0, ^| the main chaotic dynamical features of the original Lorenz 
system have been investigated through the analysis of the piecewise linear system correspond- 
ing to the choice <&(x) = \x\ and q(x) = \x\ — 1. We have thus obtained the piecewise linear 
Lorenz- like equation 

i + rjx + {\x\ - 1 + a [\x\ - lL}sgn(x) = (5) 



which has the great advantage of being exactly solvable on each side of x = 0. 
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Fig. 1 - Periodic lattice of truncated parabolae. The system evolves in the periodic constant-in-time 
potential U(x) subjected to a viscous damping and to a memory effect on the past motion. 

The Lorenz diffusion. - If in eq. (j^) one substitutes the quantities q(x) and <f>'(x) with 
x-periodic functions, the chaotic jumps between the two infinite wells of the original quartic 
potential U (x) correspond to chaotic jumps among near cells. The result is a deterministic 
diffusion in an infinite lattice, induced by a Lorenz-like chaotic dynamics. Equation ([}]) will 
be thus called the Lorenz diffusion equation. In order to use as far as possible analytical tools, 
we shall consider the unit wavelength periodic potential U(x) = ^({x} — i) (corresponding 
to q(x) — {x} — | and <&(x) = {x}), where {x} indicates the fractionary part of x. This 
potential field obviously consists of a lattice of truncated parabolae (see fig. [I]). By simple 
substitution one easily derives the equation 

x + nx + {x} - i + a[{x} =0 (6) 

where rj denotes the friction coefficient, a is the memory amplitude and (3 is related to the 
inertia whereby the system keeps memory of the past evolution. Note that rj and /3 play a 
symmetrical role in the dynamics: the solution of eq. (^) is indeed left invariant if one changes 
(3 with 77, while keeping (3(1 + a) constant. 

Inside each potential cell eq. ^ can be recasted in a third-order linear differential form. 
Indeed, by applying the operator (d/dt + (3) to each side of eq. (^), one obtains (for x =/= n) 

^3 + {[3 + r } )^ + (l+f3r ) )-+[3{l + a){x-n-^)=Q. (7) 

It is worth observing that the nonlinearity of the original model is simply reduced to a change 
of sign of the forcing term (3(l + a)({x} — \) when x crosses the cell boundaries (in our case the 
integer values). As we will see, chaotic dynamics essentially results from the unpredictability 
of the crossing times. 

Partial solutions of the third order nonlinear differential equation can be easily calculated 
inside each open interval (n, n+ 1). To obtain a global solution, such partial solutions should 
be matched at x = n by assuming that the position x, the velocity x and the memory [ |x| ]p 
are continuous, whereas the acceleration x turns out to be undefined. However, it is easily 
shown that each pair of acceleration values "immediately" before and after the crossing times 
are related by x^ +s> — x^ — sgn(x). 

The fixed points of eq. (0) are of course x = n + h , n 6 Z, and their local stability depends 
on the roots of the characteristic equation 

A 3 + (/? + ?;)A 2 + (1 + (3ri)X + /3(1 + a) = . (8) 

All the fixed points are unstable when a is larger than the critical value a c — [3~ 1 (1 + j3rf)((3 + 
77) — I. In this case eq. (||) has one real negative root (— Aq < 0) and a complex conjugate pair 
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Fig. 2 - Dependence of the exiting time from the elementary cell on the initial conditions asp, xp 
for a = 6.50, (3 = 0.19, r\ — 0.78. The graph has been obtained via numerical solution of eq. (110). 
Backward (x — n) and forward (x = n + 1) exiting cases have been distinguished by different colors. 
The structure of the graph can be guessed through the analysis of the contour lines of the function 
T = T(io,a3o) ( see re f- H] f° r further details). 



of roots with positive real part X± = A ± iu> (A > 0). For a > a c , the partial solution in the 
generic open interval (n, n + 1) can be finally written in the explicit form 

x{t) = e xt ( d cos(Lut) + C 2 sin(wf) ) + C 3 e~ A °* + n + 1/2 (9) 

where the constants C\, C%, C3 are linearly related to the (cell by cell) initial conditions. 
The motion inside each cell consists of an amplified oscillation around a central point which 
translates towards the center of the cell. A change of cell yields a discontinuous variation of 
the acceleration x and consequently of the coefficients Ci, C%, C3. 

Suppose that, at a given time, say t — 0, the particle enters the nth cell at its left 
boundary with positive velocity (the reverse case can be symmetrically analyzed). In this 
case C3 = — (Cj + 5). The question is now on whether the particle leaves the cell either from 
the left side (i.e. x — n) or from the right side of the cell (i.e. x = n+l). Assigned the model 
parameters Ao, A and o>, the minimum positive time T such that 

|e AT (CiCOs(cjT) + C 2 sm{ujT)) - {C\ + 1/2) e - A ° T | =1/2 (10) 

depends of course on C\, C%, and therefore on the initial conditions xq,xq. Unfortunately, the 
direct problem is transcendent. Moreover, as shown in fig. |2|, its solution is strongly sensitive 
to the initial conditions. This fact is a direct consequence of the crossing time definition: T is 
indeed determined by the intersection of an amplified oscillation and a decreasing exponential. 
A small change in the initial conditions may thus cause a discontinuous variation of the crossing 
time. This is the very origin of the system chaotic dynamics which suggests a stochastic 
treatment of the Lorenz diffusion equation. 

Despite the fact we have a three-dimensional space of parameters to investigate the model 
behaviors, the interesting region is actually a limited portion. This easily follows from the 
following simple considerations. For 77 large enough in eq. (^|), the motion rapidly stops in one 
of the lattice fixed points. In order to have non trivial solutions, the viscous coefficient must 
be smaller than a maximum value which can be explicitly derived from the condition a > «c 

^/(l + /? 2 ) 2 +4a/3 2 -l-/3 2 . 
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Fig. 3 - a) A typical ballistic motion with a more complex structure than a simple crossing of adjacent 
cells. Ballistic motion generally corresponds to a periodic behaviour in the cell crossing; b) trapping 
in two cells; c)-f) different kinds of trajectories generating diffusive motions. 



In the opposite limit, if r\ is too small, the friction term is negligible with respect to the 



memory term, and the resulting motion is "ballistic", i.e., ([x(t) — x(0)] 2 ) ~ t 2 (^. Analogous 
considerations can be repeated for a and f3. Inside this region of parameters the Lorenz 
diffusion equation generates a wide variety of behaviours: as we will see, the observed regimens 
are strongly sensitive to the control parameters, and furthermore this dependence is often in 
contrast with the intuitive meaning of a, (3, rj. 

In fig. U a few numerical simulation of eq. (^|) are shown, corresponding to different values 
of the parameters. These examples suggest that the variety of motion regimens ranges from 
"ballistic" ones to clearly "diffusive" , and even "trapped" in one cell or in groups of nearby 
cells. 

The analysis of the motion in the elementary cell shows that the system rapidly reaches 



( x ) Due to the deterministic nature of the system, averages should be intended over the initial conditions. 
In our numerical simulations we have typically chosen the initial conditions to be uniformly distributed over 
some real interval, and we have naturally assumed that the diffusion coefficient is independent of the choice 
of the initial ensemble. 
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Fig. 4 - a) A typical attracting set for the nth cell entering conditions xo, 'xq. The corresponding 
trajectory is shown in fig. ^e). It should be noted the presence of the straight lines defined by eq. (p"l|), 
here denoted by 7+ and 7_; b) attracting set in the plane (T, xo). The step-like structure of the set 
derives from the superposition of the points in a) on a graph like that of fig. 0. 



a steady-state. In the diffusive regimens the points corresponding to subsequent cell entering 
conditions xq , £0 are quickly attracted on a particular locus of the plane (fig. |j) . In the general 
situation this attracting set does not define a univocal map between the entering velocity and 
acceleration. However, it is not difficult to see that, for large T, xq and xq satisfy the piecewise 
linear relation || 

x - 2Xx + - (2 - A 2 - to 2 ) sgn(a-o) = . (11) 

If the diffusive regimen admits large enough times of permanence, the couple of straight lines 
defined by eq. ([ll]) therefore belongs to the entering condition attracting set (fig. [!]). Moreover, 
in the limiting case of very large Ao the attracting set exactly reduces to such lines, and in 
the steady state eq. (|ll|) is satisfied by all entering conditions Q. 

The macroscopic evolution of the system can be statistically analyzed by introducing the 
diffusion coefficient D = lim t ^oo ([x(t) — x(0)] 2 ) /(2t) . Diffusive regimens are identified by a 
finite value of D. As shown in fig. |3|, different values of the control parameters can lead to 
qualitatively very different diffusive motions. What is surprising is that the change of regimen 
is very irregular and abrupt with D, and sometimes it does not reflect our intuition on the 
role played by the parameters in eq. (|]). For instance, D shows a non-monotonic dependence 
on the viscosity rj (fig. ^), and it is possible to observe transitions from trapped (or diffusive) 
motion to ballistic ones even by increasing the viscosity itself (table |). 

Despite the somehow strange aspect of some system trajectories, it should be stressed that 
no anomalies (i.e. superdiffusive transport) have been encountered in the statistical analysis 
of diffusion, the distribution of long flies having always exponential tails. 

To conclude, we have presented a deterministic system which can generate large-scale dif- 
fusive transport induced by a Lorenz-like microscopic chaotic dynamics. Two main problems 
should still be tackled. First, the diffusion coefficient seems to have a fractal dependence on 
the control parameters, as observed for other discrete chaotic maps jl| . This fact requires more 
investigations, possibly using the analytical tools developed in ref. Second, it would be 
interesting to consider a two-dimensional extension of the Lorenz diffusion equation, mainly 
to investigate the possible emergence of anomalous diffusion induced by the increased spatial 
dimension. 
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Fig. 5 - Dependence of the diffusion coefficient D on the viscosity parameter r\ for different values of 
a and (3. 
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Table I - Different long-time behaviors as a function ofr), for fixed a and f3. Note the nontrivial, 
highly structured dependence of the diffusion coefficient on the control parameter. 
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REGIMEN 


V 


REGIMEN 


< 0.850 


ballistic 


1.425-1.725 


trapped 


0.850 - 0.875 


trapped 


1.725-1.975 


diffusive 


0.875 - 0.900 


ballistic 


1.975-2.050 


trapped 


0.900 - 1.050 


diffusive 


2.050-2.325 


diffusive 


1.050 - 1.175 


trapped 


2.325 < 


trapped 


1.175 - 1.425 


diffusive 







